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We have observed parametric generation and amplification of ultracold atom pairs. A 87 Rb Bose-Einstein 
condensate was loaded into a one-dimensional optical lattice with quasimomentum ko and spontaneously scat- 
tered into two final states with quasimomenta ki and ki . Furthermore, when a seed of atoms was first created 
with quasimomentum fci we observed parametric amplification of scattered atoms pairs in states k\ and ka, 
when the phase-matching condition was fulfilled. This process is analogous to optical parametric generation 
(OPG) and amplification (OPA) of photons and could be used to efficiently create entangled pairs of atoms. 
Furthermore, these results explain the dynamic instability of condensates in moving lattices observed in recent 
experiments. 
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Nonlinear atom optics is a novel research area born with 
the advent of Bose-Einstein condensates of alkali atoms. Un- 
like photons, ultracold atoms have a very strong nonlinearity 
directly provided by s-wave collisions, and therefore they do 
not need a nonlinear medium to provide effective interaction. 
A number of nonlinear processes first observed with photons 
have been demonstrated with matter waves such as four-wave 
mixing |1, 2], solitons |3j |4|, |5| lal, second-harmonic gener- 
ation tk. IH 191 [Toll , and sum-frequency generation 0. Non- 
linear atom optics, and in particular four-wave mixing, has 
previously been suggested as an ideal way to create entangled 
pairs of atoms ll ll flpl. However in previous four wave 
mixing experiments 1 1, 2] using condensates in free space, the 
quadratic dispersion relation for free particles only allowed for 
the phase-matching condition to be fulfilled when the magni- 
tudes of all four momenta were equal (in the center-of-mass 
frame). This is the only way in which two particles can scatter 
off each other and conserve energy and momentum. In partic- 
ular, in free space, if a condensate is moving with momentum 
ko, atoms within the condensate cannot elastically scatter into 
different momentum states, and therefore the analog to optical 
parametric generation of photons is not possible. 

The situation is very different, when an optical lattice is 
added. The lattice delivers energy in the form of the AC Stark 
effect and momentum in units of 2hkL to the atoms, where fc^ 
is the wavevector of the optical lattice. The motion of atoms in 
this periodic potential is described by a band structure, which 
deviates from the quadratic free particle dispersion curve. In a 
lattice, as recently suggested 1 13], it becomes possible for two 
atoms in the condensate to collide and scatter into a higher 
and lower quasimomentum state and conserve energy. As we 
discuss below, this can lead to dynamic instabilities of the con- 
densate, but also enables non-degenerate four-wave mixing 
and the atom-optics analog of optical parametric generation. 

Phase matching is essential for high efficiency in nonlin- 
ear processes in quantum optics including optical parametric 
generation of photons 1 14]. A modification of the dispersion 
curve has been used to demonstrated optical parametric am- 
plification in semiconductor microcavities 1 15]. In atom op- 



tics, dispersion management was used to modify the effective 
mass of atoms 1 16], and to create bright gap solitons fnll . 
Here we demonstrate that by modifying the dispersion curve 
using an optical lattice scattering processes which cannot oc- 
cur in free space become possible, and realize the matter-wave 
analogue of an Optical Parametric Generator (OPG) and an 
Optical Parametric Amplifier (OPA). 

To demonstrate the matter-wave analogue of an OPG, a 
8 'Rb Bose-Einstein condensate with quasimomentum ko was 
loaded into a one-dimensional optical lattice. To load the 
atoms at a given quasimomentum relative to the Brillouin 
zone, a one-dimensional moving optical lattice was adiabat- 
ically applied to a magnetically trapped condensate initially 
at rest in the lab frame. The lattice was created using two 
counter-propagating laser beams with frequency difference 
Sis, giving the lattice a velocity of v = -|<5^, where A is the 
wavelength of the optical lattice. In the rest frame of the lat- 
tice, the condensate has quasimomentum ko — ^Sv, where 
m is the atomic mass. By changing the detuning between the 
lattice beams, 8v, the quasimomentum ko of the condensate 
could be varied. As shown in Fig. ^1, as the value of ko 
was varied we observed elastic scattering of atom pairs into 
final states k\ and &2- The range of possible final states var- 
ied with ko due to the phase matching condition. For values 
of ko less than w .55k l the dispersion relation imposed by 
the Bloch structure of the optical lattice does not allow elas- 
tic scattering to occur. For our lattice depth of V = 0.5E rec , 
where E rec = fi 2 k\/2m, the values of ki which satisfied en- 
ergy and momentum conservation were beyond the first Bril- 
louin zone. Since the scattering process occurs within the 
first Bloch band, the atoms in state k^ have a quasimomen- 
tum ki = (ko — fci) Alod(2kT,) (see Fig. ^). As the value 
for 8v (and the resulting value of ko) is increased, the sepa- 
ration between ko and the allowed quasimomentum states fci 
and &2 decreases as clearly observed in Fig.^i. For values of 
ko above w 0.75fcz,, the final momentum states were no longer 
distinguishable, and the condensate became unstable. 

To demonstrate the matter wave analogue of optical para- 
metric amplification, we first created a small seed of atoms 
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FIG. 1: Dispersion curve for the optical lattice and experimental 
setup, (a) Band structure for a lattice depth of V=0.5E rec . The 
dashed line shows the free particle dispersion curve. The dispersion 
relation of the lattice allows two atoms with momentum fco to elasti- 
cally scatter into the final momentum states ki and ki . Energy and 
quasimomentum are conserved when ko is the average of fci and ki 
and the three points on the dispersion curve lie on a straight line. If 
ko is varied, the allowed values for fei and change. For values of 
ko below ~ 0.55fc_L, where Ul is the wavevector of the optical lat- 
tice, atoms cannot scatter elastically into different momentum states. 
The circles (squares) show allowed states ko,k\,ki for ko =0.66fc.L 
(0.70A:l). As fco is increased, the final momentum states move closer 
together. Since the scattering occurs within the lowest band of the lat- 
tice, the final momentum is ki = (2ko — fci) Mod(2fci) (b) A 87 Rb 
Bose-Einstein condensate is illuminated by two counter-propagating 
laser beams with detuning 8v, which create a moving optical lat- 
tice. The condensate is initially held at rest in a magnetic trap. In 
the rest frame of the lattice, the condensate has quasimomentum 
ko — ^-Sv. (c) As the quasimomentum fco of the condensate was 
varied, we observed elastic scattering into states fci and ki. (d) Ab- 
sorption images for different lattice detunings, 8v showing paramet- 
ric generation. After ramping up the lattice, the atoms were held for 
10 ms at a constant lattice depth. They were then released from the 
trap and imaged after 43 ms of ballistic expansion. The field of view 
is 0.5 mm x 0.3 mm 



with quasimomentum fci before ramping on the moving lat- 
tice (see Fig.|2ji. To create the seed we applied a Bragg pulse 
to the magnetically trapped condensate, outcoupling a small 
fraction of atoms into the momentum state ksragg 11811 . Im- 
mediately after applying the pulse, the optical lattice was adi- 
abatically ramped on. In the rest frame of the lattice, the seed 
has quasimomentum fci = ksragg + ko- The phase matching 
condition for a given seed momentum fc Bragg was found by 
varying the frequency difference 5v of the lattice, and there- 
fore the quasimomenta fco and fci of the atoms. As shown in 
Fig. |2jl, when the phase matching condition was fulfilled, we 
observed amplification of the seed fci as well as its conjugate 
momenta k^. The growth of the seed and its conjugate mo- 
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FIG. 2: Parametric amplification of scattered atom pairs in a ID 
optical lattice, (a) First, a 2 ms Bragg pulse was applied to the con- 
densate, (b) The Bragg pulse seeded atoms along the long axis of 
the condensate with momentum ksragg = (k a — fci,) in the lab frame 
(c) The optical lattice was then adiabatically ramped on and applied 
for 10 ms. When the phase matching condition was fulfilled, para- 
metric amplification of atoms in the seeded state fci and its conjugate 
momentum state ki was observed, (d) Resonance curve showing 
amplification of ki, when fci is seeded. Amplification occurred only 
when the phase-matching condition was met. For a fixed fc Bragg, 
the resonance condition was found by varying the detuning 8u of the 
lattice. The data was taken for ksragg = 0.43fc,r,. The fraction of 
amplified atoms was obtained by subtracting images with and with- 
out the seed pulse, (e) Absorption images showing amplification of 
fci and ki when the phase matching condition is met. The center of 
the resonance was at 8v ~ 5450Hz, close to the calculated value of 
5v » 5350Hz. 



mentum as a function of time are shown in Fig. [3] 

The experiments were performed using an elongated 87 Rb 
Bose-Einstein condensate created in a cloverleaf-type Ioffe- 
Pritchard magnetic trap previously described in Ref. 11911 . 
The magnetic trap had a radial (axial) trap frequency of 35(8) 
Hz. The condensate, containing between 0.5 — 3.0 x 10 s 
atoms, was produced in the \5?Su2,F — l,mj? = — 1} 
state. The Bragg pulse was created with two laser beams de- 
rived from the same laser, which was red-detuned from the 
5 2 Si /2 , F = 1 -> 5 2 P 3 / 2 , F = 1 transition at A = 780nmby 
400MHz, and was 7r-polarized. As shown in Fig |2j the Bragg 
beams were aligned such that atoms were outcoupled along 
the long axis of the condensate. The intensity of the Bragg 
pulse was chosen such that less than 5% of the initial conden- 
sate was outcoupled into the state fc Bragg, and the length of 
the pulse was 2ms. The angle between the Bragg beams could 
be varied to change the momentum of the outcoupled atoms. 
The optical lattice was created using two counter-propagating 
beams derived from the same laser with A = 1064nm, and the 
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FIG. 3: Growth curve for atomic population in quasimomentum 
states k\ and &2 when the process is seeded, (a) Amplification of 
atoms with quasimomentum fci (solid points), and with the conju- 
gate momentum ki (open points), when state fei was seeded, (b) 
Amplification of atoms with quasimomentum ki (solid points), and 
with the conjugate momentum (open points), when state &2 was 
seeded. The values for k^,k\,k% were 0.66fcL,0.23/ct,and -0.92/cl 
respectively. The solid lines shows the expected gain using Eq. Q 
with a variable scale factor as the only free parameter. 



frequency of the two beams were controlled by two separate 
acousto-optic Modulators (AOMs) driven with a frequency 
difference 5v. The lattice was also aligned along the long axis 
of the condensate, and was ramped on in 1 ms using an expo- 
nential ramp. After the condensate was held in the lattice for a 
variable time r it was then released from the trap and imaged 
after 43 ms of ballistic expansion. 

For all of our experiments, the depth of the optical lattice 
was V = 0.5 E rec with a band structure shown in Fig. [2 
When the process was not seeded, atoms were elastically scat- 
tered into a narrow band of states k\ and k%, where both en- 
ergy and momentum was conserved. However the popula- 
tion in neither state was large enough for amplification to be 
observed. When the process was seeded, amplification oc- 
curred when the quasimomentum was tuned such that energy 
and momentum were conserved for the states ko, fci, and the 
conjugate momentum k 2 . In our experiment, the difference 
Afc = ko — ki between the quasimomentum of the condensate 
fco and seed fci was set by the angle of the initial Bragg pulse. 
In the rest frame of the lattice, for a given frequency differ- 



ence 8v the quasimomentum of the atoms was ko — ^^5v, 
and that of the seed was k\ = k Bragg + ko. 

For a given Bragg angle, there is only one set of quasimo- 
menta ko,ki, and k 2 where the phase-matching condition is 
fulfilled. To find this point, we varied the velocity of the mov- 
ing lattice for fixed hold times. Results for fcs ragg =0.43fcL are 
shown in Fig. [2] The phase-matched value for k 2 is at 1.08/cl, 
beyond the boundary of the first Briollouin zone. Therefore, 
the atoms are observed with a momentum k2 = -0.92fc^. For 
ksragg = 0.43/cL,0.34fc£ and 0.28/cl, we observed resonances 
at (5i/=5450Hz, 5750 Hz and 6100Hz respectively. For these 
Bragg angles and our lattice depth, the expected values were 
5350Hz, 5700Hz and 6050Hz. 

In Fig |3j 5% of the initial condensate containing No = 
1.3(3) x 10 5 atoms was outcoupled with ksragg = 0.43fc£. The 
gain for the process is determined by the strength of the non- 
linear interaction U between atoms in the condensate. For a 
shallow lattice U = - , where a is the s-wave scattering 
length. From this we can estimate the maximum amplifica- 
tion rate to be r\ = 2nolI/h H, with Ni(2) = T]N 2 (i), where 
Ni( 2 ) is the number of atoms in the momentum state ku2)- 
For Nq = 1.3(3) xlO 5 , the maximum growth rate should be rj 
= 540 Hz. The amplification rate will decrease as the state ko 
is depleted. However, for our small seeds, the amplification 
was limited by the loss of overlap between the condensate and 
the amplified pulses. The Thomas-Fermi radius (Rtf )of the 
condensate in the axial direction was 33/im, and the recoil ve- 
locity (tVec) for the final states k\ and k 2 with respect to the 
initial condensate was u rec =1.8 /im/ms and 6.8/im/ms respec- 
tively. The overlap integral between the amplified atoms and 
the initial condensate can be approximated as a Gaussian with 
time constant r c w 0.75RTF/v re c, which for our parameters 
is 3.75 ms. We compare our results to the modified rate equa- 
tion 

-t 2 

#3(1) =^i(2) e^" (1) 

Since atoms are scattered into states k\ and k 2 in pairs, one 
would expect that the final atom number in the two states (mi- 
nus the initial seed) are equal. Instead, we observe a smaller 
number in state k 2 which we ascribe to the proximity of k 2 
to the boundary of Brillouin zone. This leads to instabilities, 
where atoms in state k 2 are scattered into other momentum 
states or into higher bands. If we allow a variable scale factor 
to correct for this in our model, as shown in Fig. [3] the final 
values for N\^ are in agreement with the experimental data. 

Amplification was also observed when atoms were seeded 
in state k 2 . Due to the geometry of our experimental setup, we 
were unable to load atoms directly into k 2 = -0.92fc/,. How- 
ever when atoms with quasi-momenta k = 1 .08fcL were loaded 
into the lattice, the ramp-up was no longer adiabatic due to 
their proximity to the boundary of the first Brillouin zone. Be- 
cause of this, atoms from the seed were loaded into both the 
second Bloch band (with k = 1.08fcL) and the ground state 
(with k = -0.92fci = k 2 ). As shown in Fig. the gain for 
this process was almost identical to when atoms were seeded 
in state ki . 
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The loss of overlap could be alleviated by using a more 
extreme trap geometry in which the condensate is more elon- 
gated, e.g. by confining atoms in a tight transverse optical 
lattice. In this configuration, it may be possible to observe the 
parametric scattering dynamics for longer time scales, which 
may allow for the observation of Rabi oscillations between 
fco and k\, hi as predicted in Ref. 11 311 . For longer coher- 
ence times, parametric amplification could also be an efficient 
means of producing pairs of momentum entangled atoms for 
quantum information applications 11 UM 211 . 

For high atom numbers, and for large values of fco the 
condensate became unstable, and scattered into a broad band 
of final momentum states (Fig. |2^). For ksragg = 0.43/cl, 
the energy of atoms outcoupled by the Bragg beams was 
w 370H z, whereas the chemical potential of the condensate 
was ss 300-ffz. Because of this, if the atom number was in- 
creased significantly the momentum peaks were no longer dis- 
tinguishable. When the chemical potential of the condensate 
was larger than the separation between the phase matched mo- 
mentum states, the process was self-seeded, i.e the momentum 
spread of the initial condensate contained atoms with momen- 
tum fco,fci,and &2, and considerable scattering occurred. Sim- 
ilarly, if the atom number was kept constant, and the value of 
fco was increased, the phase-matched momentum states move 
closer together until they are no longer distinguishable. This 
occurred at values of fco above w 0.75fci, and we observed 
a dynamic instability. For larger atom numbers, the critical 
value of fco decreases. For values of fco less than 0.55fco elas- 
tic scattering cannot occur, and the system should be stable 
for all atom numbers. Dynamic instabilities of condensates 
in moving lattices were recently observed in Ref. 1 20, 21 ]. In 
Ref. 112 ill , the chemical potential was a factor of 3 higher than 
in our experiment, leading to a dynamic instability for all val- 
ues of fco above 0.55fcL- Although discrete momentum states 
could not be observed in those experiments, it is likely that 
the mechanism for the dynamic instability is self-seeded para- 
metric amplification. 

In conclusion, we have demonstrated a matter wave ana- 
logue of both optical parametric generation and optical para- 
metric amplification using a condensate moving in a one- 
dimensional optical lattice. The optical lattice modified the 
dispersion curve and ensured phase matching. If the separa- 
tion of the phase matched momentum states becomes less than 
the speed of sound, a condensate will self-seed the process and 
become dynamically unstable. 

The authors thank K. M0lmer for pointing out this physi- 
cal process to us, M. Inguscio for insightful discussions, M. 



Saba for insightful discussions and a critical reading of the 
manuscript, and R Medley for experimental assistance. This 
work was supported by NSF and ARO. 



* URL: |http : / / cua .mit . edu/ketterle_group/| 
[1] L. Deng et al, Nature 398, 218 (1999). 

[2] J. M. Vogels, K. Xu, and W. Ketterle, Phys. Rev. Lett. 89, 

020401 (2002). 
[3] J. Denschlag et al., Science 287, 97 (2000). 
[4] S. Burger, K. Bongs, S. Dettmer, W. Ertmer, K. Sengstock, 

A. Sanpera, G. V Shlyapnikov, and M. Lewenstein, Phys. Rev. 

Lett 83, 5198 (1999). 
[5] B. P. Anderson et al., Phys. Rev. Lett 86, 2926 (2001). 
[6] K. E. Strecker, G. B. Partridge, A. G. Truscott, and R. G. Hulet, 

Nature (London) 417, 150 (2002). 
[7] R. Wynar, R. S. Freeland, D. J. Han, C. Ryu, and D. J. Heinzen, 

Science 287, 1016 (2000). 
[8] J. R. Abo-Shaeer, D. E. Miller, J. K. Chin, K. Xu, 

T. Mukaiyama, and W. Ketterle, Phys. Rev. Lett 94, 040405 

(2005). 

[9] P. D. Lett, K. Helmerson, W. D. Phillips, L. P. Ratliff, S. L. 
Rolston, and M. E. Wagshul, Phys. Rev. Lett 71, 2200 (1993). 
[10] J. D. Miller, R. A. Cline, and D. J. Heinzen, Phys. Rev. Lett 71, 
2204 (1993). 

[11] L.-M. Duan, A. S0rensen, J. I. Cirac, and P. Zoller, Phys. Rev. 

Lett 85, 3991 (2000). 
[12] H. Pu and P. Meystre, Phys. Rev. Lett. 85, 3987 (2000). 
[13] K. M. Hilligs0e and K. M0lmer, Phys. Rev. A 71, 041602(R) 

(2005). 

[14] E. Rosencher, Optoelectronics (Cambridge University Press, 

Cambridge, UK, 2002). 
[15] P. G. Savvidis, J. J. Baumberg, R. M. Stevenson, M. S. Skol- 

nick, D. M. Whittaker, and J. S. Roberts, Phys. Rev. Lett. 84, 

1547 (2000). 

[16] B. Eiermann, P. Treutlein, T. Anker, M. Albiez, M. Taglieber, 
K.-P. Marzlin, and M. K. Oberthaler, Phys. Rev. Lett 91, 060402 

(2003) . 

[17] B. Eiermann, T. Anker, M. Albiez, M. Taglieber, P. Treutlein, 
K.-P. Marzlin, and M. K. Oberthaler, Phys. Rev. Lett 92, 230401 

(2004) . 

[18] J. Stenger, S. Inouye, A. P. Chikkatur, D. M. Stamper-Kurn, 
D. E. Pritchard, and W. Ketterle, Phys. Rev. Lett. 82, 4569 
(1999). 

[19] D. Schneble, Y. Torii, M. Boyd, E. Streed, D. E. Pritchard, and 

W. Ketterle, Science 300, 475 (2003). 
[20] L. Fallani, L. De Sarlo, J. E. Lye, M. Modugno, R. Saers, 

C. Fort, and M. Inguscio, Phys. Rev. Lett 93, 140406 (2004). 
[21] L. De Sarlo, L. Fallani, J. E. Lye, M. Modugno, R. Saers, 

C. Fort, and M. Inguscio, Phys. Rev. A 72, 013603 (2005). 



